Polarization-Entangled Light Pulses of 10^ Photons 



(N 

O 

(N 

OO 



> 

o 

<N 



X: 



Timur Sh. Iskhakov/ Ivan N. Agafonov,^ Maria V. Chckhova/'^ and Gcrd Leuchs^-^ 

^ Max Planck Institute for the Science of Light, Gunther-Scharowsky-Strafie 1/Bau 24, 91058 Erlangen, Germany 
^Department of Physics, M.V.Lomonosov Moscow State University, 
Leninskie Gory, 119991 Moscow, Russia 
^University of Erlangen-Niirnberg, Staudtstrasse 7/B2, 91058 Erlangen, Germany 

We experimentally demonstrate polarization entanglement for squeezed vacuum pulses containing 
more than 10^ photons. We also study photon-number entanglement by calculating the Schmidt 
number and measuring its operational counterpart. Theoretically, our pulses are the more entangled 
the brighter they are. This promises important applications in quantum technologies, especially 
photonic quantum gates and quantum memories. 
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Entanglement is the signature of the quantum world. 
One part of an entangled system has its properties fully 
undefined yet fully correlated with the properties of 
its counterpart Can this behavior be observed for 
large objects? Recently, entanglement was discovered for 
macroscopic material systems 0, It is tempting to 
observe it for bright photonic states because bright 
light is much more efficient in interactions than single 
photons. For bright squeezed vacuum, very different from 
usual squeezed light, entanglement was discussed theo- 
retically [H but never tested experimentally. Coin- 
cidence measurements could only reveal entanglement for 
up to 12 photons [ll 



The clue to the observation of entanglement for bright 
squeezed vacuum is in re gist ering, instead of single pho- 
tons and coincidences [ol. Illl [l2| , fluctuations of macro- 
scopic intensities and measuring the variances of intensity 
differences A great advantage of this measurement 
is that it is robust against the multimode detection. In 
our experiment, by applying this technique to entangled 
bright squeezed vacuum [14| we test a multimode sepa- 
rability condition formulated in terms of variances. 

The states of entangled squeezed vacuum, also known 
as macroscopic Bell states, can be obtained via para- 
metric down-conversion (PDC) in two type-I nonlinear 
crystals For example, the singlet state is generated 
by the Hamiltonian [1 i, i, [ll| 

(1) 



where a\ arc photon creation operators in beams A, B, 
which in our case have the same direction but different 
wavelengths Xa,^b- The subscripts H,V stand for the 
horizontal and vertical polarization, and the parameter 
G depends on the crystal properties and the pump power. 
The state can be written as a Fock-state expansion [l[ , 
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Clearly, photon numbers in beams A and B are ex- 
actly the same. Moreover, if polarization beamsplitters 
are placed into the beams, the number of transmitted 
photons in beam A will be exactly equal to the number 
of reflected photons in beam B, and vice versa (Fig. [T|). 
Correlations will be maintained at any orientations of 
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FIG. 1. Quantum correlations of the macroscopic singlet Bell 
state. The source S emits light pulses into beams A and B 
(spatially separated for clarity). In each pulse, photon num- 
bers emitted into any orthogonal polarization modes in the 
two beams are random but exactly equal. 



the polarizers, or quarter-wave plates inserted in front of 
them, as long as they are the same for both beams. This 
is because the operator expression in I'l'^') is invari- 
ant to any polarization transformation. Note that such 
perfect correlations are only manifested by squeezed vac- 
uum and not by displaced squeezed states, which contain 
a huge coherent component and only a small proportion 
of squeezed vacuum. 

Such polarization correlations are similar to the ones 
manifested by two-photon Bell states but involve far 
larger photon numbers. Besides their fundamental in- 
terest, they are important for quantum information pro- 
tocols based on light-light (quantum gates) and light- 
matter (quantum memory) interactions. One can men- 
tion quantum memory based on bright squeezed vac- 
uum [15| and up-conversion of such states [l6| . The latter 
suggests a new field of research, nonlinear optics of en- 
tangled states. 

Separability condition. How can we prove in ex- 
periment that the state PH)) is entangled? According to 
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Ref . [H , if a bipartite system containing two macroscopic 
light beams A, B (Fig. [1]) is separable, it satisfies a cer- 
tain condition. Violation of this condition indicates that 
the state is non-separable (entangled if it is pure) . 

To simplify comparison with experiment, we derive a 
necessary condition of separability in terms of the Stokes 
parameters and their variances |17| . This approach en- 
ables us to prove a stronger condition than the one of 
Ref. [ij. It is important that our consideration is also 
valid for multimode beams. 

As shown in Section A of the supplementary material, 
for a separable state, the sum of the three Stokes vari- 
ances AS'j^. i = 1,2,3, cannot be smaller than twice the 
total photon number (S'o) , 



^ASf/(^o)>2. 
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A similar inequality was proved for atomic ensembles in 
Ref. 0. 

Inequality ([3]) is often mentioned as one of the un- 
certainty relations in polarization quantum optics (see, 
for instance, 0, [l^). Indeed, it follows directly from 
the well-known equality for the Stokes operators Si Q, 
Si + Sl + Si ^ 5*0 (5*0 + 2). It should be noted, how- 
ever, that this operator equality holds true only in the 
case where, apart of the two polarization modes, the 
light beam contains only a single frequency and angu- 
lar mode [21, 22 1 . Thus, inequality ([3|) is not of general 
validity. In fact, it is a necessary condition of separabil- 
ity. Its violation indicates that a beam is non-separable, 
i.e., is a sufficient condition of non-separability. As we 
show below, Eq. ([3]) is violated in our experiment. 

The experiment was performed with the macro- 
scopic singlet Bell state \^mac)^ similar to the one con- 
sidered in [llUJS, 11 1 . The setup (Fig. is described in 
detail in Refs. ja[l7|. 

Theoretically, the singlet state | '^mlc) has three Stokes 
parameters equal to zero, ( 51.2, 3) = 0, as well as the 
correspon ding variances, AS^ 2 3 — Oj higher-order 
moments [2l|. Thus, in theory condition (U) is always 
violated, as its left-hand side is zero. In practice, achiev- 
ing a zero variance of any Stokes observable is impossible. 
The noise is caused by the inevitable losses (including the 
non-ideal quantum efficiency of the detectors) and the 
imperfect mode matching. To optimize the mode match- 
ing, beams A, B are filtered in the angle separately. 

Testing condition ^ requires the measurement of vari- 
ances for three Stokes observables and the total pho- 
ton number {So), which is the shot-noise level. The 
variances of >S'i_2,3 and the mean photon number (S'o) 
were measured by analyzing the statistics over 20000 
pulses [2^. Typical photon numbers per pulse were 
10^ (Fig. [2]). This was due to a very large number of 
modes collected. It is known that collincar type-I phase 



FIG. 2. Top: the experimental setup. Orthogonally po- 
larized squeezed vacuums are generated in two type-I BBO 
crystals and overlapped at a polarizing beamsplitter (PBS); 
the residual pump is eliminated by dichroic mirrors (DM) and 
a long-pass filter OG. The interferometer is balanced using a 
trombone prism. A dichroic plate (DP) is inserted for pro- 
ducing the macroscopic singlet state. Two apertures (A1,A2) 
placed in the focal plane of a lens select the angular spectra of 
the beams at two wavelengths, separated by a dichroic mirror 
DM1 and joined together by another dichroic mirror DM2. 
The measurement part also includes a Glan prism (GP), a 
half- wave plate (HWP) or a quarter- wave plate (QWP), and 
two detectors. Bottom left: number of photons per pulse ver- 
sus the pump power. Bottom right: variance of the Stokes 
observable plotted versus the direction in 3D (the object in- 
side the sphere). The outer sphere shows the shot-noise level. 



matching is characterized by a large number of angu- 
lar Schmidt modes (23 |. By accepting, with our angular 
apertures, nearly whole angular spectra at wavelengths 
AajAb, we collected about 10^ angular modes and 10^ 
frequency modes Q. At the same time, the number of 
photons per mode was mesoscopic. The bottom-left part 
of Fig. [2] shows the output /input characteristic of the 
down-converter; the dependence is almost linear and its 
fit yields the maximum gain 0.33 ±0.06 corresponding to 
the number of photons per mode 0.12 ± 0.04. However, 
condition ([3]) is invariant to the number of modes [l7| . 
This is why it is suitable for testing multimode states; on 
the contrary, traditional Wigner-function measurement 
requires single-mode states and is therefore not applica- 
ble here. Besides, measurement of the photon-number 
variance for squeezed vacuum has been shown to be in- 
variant to the gain, at least up to values F ~ 2 

Figure [3] shows the left-hand side of inequality ([3]) plot- 
ted against the diameter of the Al aperture, Di. For all 
points below the dashed line, the necessary condition of 
separability is violated. We see that with the transverse 
modes properly matched, it is violated by more than 5 
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standard deviations. 




FIG. 3. The left-hand part of the separability condition ([3|) 
versus the diameter of aperture Al, the other aperture diam- 
eter being 8.9 mm. As expected, the best noise suppression is 
observed for Di satisfying the condition D1/D2 = Xa/Xb 0, 
for which the mode matching is optimal. The dashed line is 
the boundary set by the separability condition. 



Hence, we have experimentally demonstrated that 
the macroscopic singlet Bell state is polarization non- 
separable. Note that it is prepared pure, and only im- 
perfections of the detection setup make it mix with the 
vacuum. A pure non-separable state should be able to 
violate Bell's inequalities but apparently a new form 
of Bell's inequality should be derived for this state. 

Photon-number entanglement of the singlet Bell 
state (fTOj) can be characterized by noticing that the state 
can be rewritten as a product of two Schmidt decompo- 
sitions in the Fock basis. 
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Here, A„ 

means a Fock state in beam A with n photons in the hor- 
izontally polarized mode. The notation for beam B and 
mode V is similar. Clearly, the state can be represented 
as a product of two entangled states, one of them involv- 
ing modes AH, BV and the other one, modes AV, BH . 

For each of the states |vl/i.2), the Schmidt number is 
Ki = K2 = 1 + 2iVo, where Nq = sinh^ F is the mean 
photon number. The total Schmidt number is their prod- 
uct. Note that none of the states |*I'i,2), taken sepa- 
rately, violates condition ([3]), although they both mani- 
fest photon- number entanglement. This shows that there 
is a difference between polarization and photon-number 
entanglement of macroscopic Bell states. 



If there are many (M) independent frcquency- 
wavevector mode pairs, each containing a state of the 
form ([T0| , then the total Schmidt number is given by the 
product, K = (l-|-2A'o)^^'^, which is extremely large. Un- 
der certain experimental conditions this huge amount of 
entanglement could be used. However, in our experiment 
we treat the whole ensemble of modes jointly; moreover, 
the detection scheme also combines the two wavelengths. 
With the only partition being the polarization one, the 
Schmidt decomposition for our state can be written as 
for a two-mode squeezed vacuum, but with A„ given 
by the Poissonian distribution (soj with the mean value 
N = 2MNq. Then the Schmidt number can be calculated 



as K 



//o(27V), where Iq is the zero-order modified 



Bessel function of the first kind. At large N, K w 2 v t:N . 
In our experiment, this yields K Ki 10^. This is much less 
than in the single-mode case because we do not address 
each mode separately and only deal with their ensemble. 
We stress that the Schmidt number is not an operational 
measure as it cannot be directly measured in experiment. 
However, there exists its operational counterpart. 

Operational measure. The Schmidt number can 
be understood as the effective number of the Schmidt 
modes. For continuous variables like wavevectors or fre- 
quencies of single photons, its operational counterpart 
was proposed [26| as the ratio of the unconditional width 
of one subsystem w.r.t. some parameter to the width of 
the corresponding conditional distribution. 

By analogy, consider the following procedure. In ex- 
periment, we obtain the photon- number probability dis- 
tribution in, say, beam B {unconditional distribution). 
Next, we measure the photon- number conditional distri- 
bution for beam B by post-selecting only those pulses for 
which the photon number in beam A is fixed, or within 
certain narrow bounds. The ratio of the widths for the 
two distributions can be considered as an operational 
measure of entanglement. Note that the narrowing of 
the photon-number conditional distribution is a typical 
feature of twin beams [2^ but has never been considered 
as a measure of entanglement. 

For each pair of correlated modes (for instance, AH — 
BV), the joint probability distribution contains a factor 
SnA,nB (Fig- UK)- Therefore, while the unconditional dis- 
tribution for beam B has a negative exponential shape, 
the conditional probability is only nonzero for a sin- 
gle value of n^. In the multimode case, the uncon- 
ditional distribution is Poissonian (Fig. HJj) while the 
conditional distribution is again of unity width. This 
gives R = 2y/2N\n2, almost coinciding with the Schmidt 
number. Unfortunately, the conditional distribution gets 
broadened due to losses, and the resulting R value be- 
comes Rrf — 1/^/1 — 77, where 77 is the overall detection 
efficiency. Thus, the accessible degree of entanglement is 
given by only the detection efficiency and turns out to be 
much less than the Schmidt number. 

In our experiment, the unconditional distribution is 
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broadened due to the excess noise of the pump. Therefore 
we find the numerator of the R ratio by measuring the 
width of a shot-noise hmited source with the same mean 
photon number (Fig. Hh, red sohd hne). The degree of 
entanglement is then measured to be 1.53 ± 0.05. This 
corresponds to the detection efficiency rj = 0.57. 




c 






—^^Shot noise limited . 




/\ 




// A Conditional 




/' \ 




/' A 




/ , R=1.53 \ 
/ ^ \ 







FIG. 4. Photon-number distributions, (a) single-mode case: 
the joint probability distribution (blue line in the plane ua ~ 
ub), the unconditional distribution (exponentially decaying 
red line), and the conditional distribution (green dashed line), 
(b) the same distributions in the case of a realistic experi- 
ment (many modes, non-ideal detection), (c) experimentally 
obtained distributions: the unconditional one, with the ex- 
cess noise eliminated (red solid line) and the conditional one 
(green dashed line). The electronic noise is subtracted. 



The triplet states. We have proved the violation of 
the separability condition ([3]) for the macroscopic singlet 
Bell state. Since the triplet states can be obtained from 
the singlet one via local unitary transformations, they 
are entangled as well. The separability conditions for 
these states can be derived from condition ([3]) via the 
corresponding unitary transformations [29 ] . 

The Schmidt decomposition for the triplet states is 
similar to the one for the singlet state, and the Schmidt 
number is the same. The only difference is that for the 
singlet state, the polarization modes can be chosen in any 
way. For each of the triplet states, there is a unique choice 
of polarization modes to see entanglement: it should be 
horizontally and vertically polarized modes for |^'^,|ic), 
diagonally polarized modes for |$Lac), and right- and 
left-circularly polarized modes for \^m}c) jsj. 

In conclusion, we have tested the macroscopic singlet 
Bell state, containing two beams of different wavelengths, 
for separability. The results convincingly prove that the 
state is non-separable with respect to polarization ob- 



servables. As the photon numbers per pulse are as high 
as 10^, and the state is prepared pure, this can be con- 
sidered as a proof of macroscopic entanglement. As a 
measure of photon-number entanglement, we have calcu- 
lated the Schmidt number, which turned out to be given 
by the total number of photons in a pulse. Theoretically, 
our multiphoton state is highly entangled. The entan- 
glement can be confirmed by the measurement of the 
photon-number distribution for one beam (unconditional 
distribution) and the distribution conditioned on regis- 
tering a certain number of photons in the other beam 
(conditional distribution) . As the former is broader than 
the latter, we conclude that the state is photon-number 
entangled. The measured degree of entanglement and 
separability condition violation are both reduced due to 
the inefficient detection. However, there is apparently a 
difference between photon-number and polarization en- 
tanglement. 
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Supplementary information: 
Polarization-Entangled Light Pulses of 10^ 
Photons 



Derivation of the separability condition. 

Li this derivation, we will proceed along the same line 
as Simon and Bouwmeester with the difference that 
we will use polarization variables. This makes the proof 
more consistent with the experimental part. Another im- 
portant difference is that we go further than the authors 
of Ref. 1] and prove a stronger inequality, bounding from 
below not the mean square value of a vector operator 
but its variance. Finally, we extend the inequality to the 
multi-mode case. 

Consider a bipartite system containing two macro- 
scopic light beams, A and B, which are distinguishable in 
wavelength (our case) or wavevector (the case considered 
in Ref. 1]). The Stokes vector of such a system, de- 
fined as a set of three Stokes operators, S = {^i; S'2; 5*3} 
is given by the sum of the Stokes vectors for beams 
A,B [l|,[3, 



S = 



(4) 



Suppose that the beams are separable. Then the den- 
sity matrix of the total state is p = Pjpf pf ■ The 

'3' 



square of the Stokes vector operator is = Sf + S^ + S^, 



and its mean value can be written as 
(S2)=5]p,([S-4]2),+5]p,([S^]2),+2 5^p,(S-4),(S^ 



The subscript j by the angular brackets means that the 
averaging is made over the state pf'^- Here, pj > are 
classical probabilities, ^jPj = 1- It is important that 
in the last term, the averaging goes separately for each 
subsystem due to separability, (S'^S'^)j = (S^)j(S^)j. 
The variance of S takes the form 

as' = (S2) - ((S^ + S^))2 = 

^p,([S^]2),+^p,([S^]2),+2^p,(S-4)^(S%- 

j j 



The last three terms here are squared mean values cal- 
culated by averaging with usual classical probabilities 



6 



Pj > 0. This implies the positivity of the variance, and 
for each term we can write 



<Y.P^(^k + Si)l A; = 1,2,3, 



which leads us to an inequality for the variance. 



AS > 



i 



Ai2\ 



j 



Bi2\ 



2E 



i 



pAs^ + s?) 



B\2 

j 



3 

Again, due to separability, {SfSi)j + /j 



In the special case of light unpolarized in the first order, 
for which (S) = 0, both conditions coincide. 

It should be emphasized that inequality ©, as well 
as inequality ([7]), follows from the separability assump- 
tion and is therefore a necessary condition of separability. 
Hence, the opposite inequalities are sufficient conditions 
of non- separability and can be called non- separability wit- 
nesses. 

In the multimode case, where each beam A, B contains 
a set of independent modes, one can write inequality ^ 



inequality 



i 

This leads to the 



cS\2for each pair of modes Ak,Bk, 
^3 / j- 

ASl > 2{So)k. 



AS > 



E 



Ai2\ 



i 



B-\2\ 



E 



To each beam A, B taken separately, we can apply the 
well-known operator equality Q , 



1A.S12 



(5) 



which yields 



AS^ > 



-2)),+Y,p,{S^{S^ + 2)), 
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Note that up to this step, no assumption was made 
that the beams A, B are single-mode. However, equality 
([S]) is only valid for a single-mode beam. Hence, so far 
we restrict our consideration to a single pair of modes. 
The multi-mode case will be considered below. 
Since 



(S^'^)2 



3 \-U O 

is the first-order degree of polarization 
can be rewritten as 



where < Pf'^ < 1 



the inequality 



as' > 5:^,00^(^0^ + 2)), 



E^ 



pAS^{S^+2)), 



^P3{s^)',-Y.pAS^)1 



3 



Note that {[S^A^) > (M) 



negative), and J^jPji^o + S^)j 
the separability condition in the form 

AS2 > 2(5o). 



(the variance is non- 
= (So), hence we get 



(6) 



A similar inequality was proved for separable atomic en- 
sembles in Ref. Q. This inequality is stronger than the 
separability condition derived in Ref. 1], 



(S^) > 2{So) 



(7) 



(8) 



The variance of any observable for an ensemble of inde- 
pendent modes is equal to the sum of separate variances, 
AS^ = X^fc^^L ^^"^ same is true for the mean val- 
ues, (5*0) — J2ki^'^)k- Then, by summing inequalities 
we again obtain inequality ^ but this time it is valid for 
a multi-mode system. Therefore, the opposite inequality. 



AS2 < 2(5*0) 



(9) 



is a witness of non-separability for any kind of states, 
including mixed and multi-mode ones. 



Preparation of the I"!* 



{-) 



state 



We pump two orthogonally oriented type-I BBO crys- 
tals placed in two arms of a Mach-Zehnder interferom- 
eter (MZI) with a strong pulsed beam (third harmonic 
of Nd:YAG laser with the pulse duration 18 ps, pulse 
energy up to 0.2 mJ, and repetition rate 1 kHz). The 
four- mode state is realized by, first, using frequency non- 
degenerate phase matching, signal and idler wavelengths 
being — 635 nm and Xb — 805 nm, respectively, and, 
second, overlapping two orthogonally polarized squeezed 
vacuum beams on a polarizing beamsplitter. The phase 
between the beams is set to be tt using a piezoelectric 
feed; it is important that the MZI is balanced with an 
accuracy better than the coherence time of the pump, 
which is 5 ps. To provide this, we put one of the inter- 
ferometer mirrors on a piezoelectric feed combined with 
a micrometric translation stage. Fig. [5] shows how the 
interference visibility depends on the 'rough' shift of the 
micrometric stage. The interference is observed by mea- 
suring the 5*2 Stokes variance versus the bias voltage ap- 
plied to the piezoelectric feed. The bias voltage changes 
the phase between the squeezed vacuums produced in 
the MZ interferometer, so that the state at the output is 

|^'((p)) = exp[r(a]j6|^ + e^'^albl) + h.c.]|vac). (10) 



At phase ip — n, the state us |$Lac) 
fiuctuations of the 5*2 observable [2| . 



, and it has reduced 
At iy9 = 0, the state 
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FIG. 5. Interference pattern envelope (minimum and max- 
imum NRF values for the S2 Stokes observable) versus the 
path length difference between the arms of the MZ interfer- 
ometer. 
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ing a separate iris aperture for each wavelength. Both 
apertures are placed in the focal plane of a lens with the 
focal length 300 mm. The crystals, in their turn, are 
placed at a distance of 300 mm from the lens. This way, 
the aperture diameters can be adjusted independently so 
that the same range of transverse wavevector (the same 
numbers of transverse modes) for both signal and idler 
wavelengths can be filtered out. The 635 nm and 805 
nm beams are then joined on another dichroic mirror, 
and for the whole beam the Stokes measurement is per- 
formed. Note that when two beams of different colors 
are registered by the same detectors, the Stokes observ- 
ables simply add up, as the detectors do not resolve the 
frequency beats. 

Our detectors are p-i-n diodes followed by charge- 
sensitive amplifiers. For each light pulse, they both pro- 
duce output pulses with a fixed shape and the amplitude 
scaling linearly with the total number of registered elec- 
trons. The output pulses are time-integrated by means 
of an AD card. The electronic noise is mainly caused by 
the electric circuit as well as the digitization noise of the 
AD card and amounts to 300 electrons/pulse. Before the 
measurement, the detectors are calibrated using a shot- 
noise limited source (strongly attenuated laser radiation), 
by measuring the variance of the difference signal S- as 
a function of the mean sum signal (5+). In the range 
of (5+) on the order of 10^ - 10^ photons per pulse, the 
dependence is linear, AS'i = a(S'+) -I- A5'g, A5g being 
the electronic noise contribution. The first term of this 
dependence is considered as the shot-noise level and used 
in the calculation of NRF. 



FIG. 6. NRF of the Stokes observable S2 as a function of 
the bias voltage at the piezoelectric feed. The arrow shows 
the point where the fluctuations of S-2 are suppressed, which 
indicates the production of |<I>kac) and hence, the phase = 

TT. 



is |$^2c) and the fluctuations of 5*2 are enhanced. Scan- 
ning of the bias voltage and hence the phase <~p creates a 
sine variation of the normalized variance (noise reduction 
factor, NRF) of S2 (Fig.[6|). By plotting the minimal and 
maximal NRF values versus the reading of the micromet- 
ric stage we get the envelope of the interference pattern 
(Fig.©. 

The residual pump is eliminated using a dichroic mir- 
ror in each arm and a colour-glass filter after the MZI. 
The output state |$m2c) becomes l^mic) in the 45°-tilted 
basis. It is turned into the singlet state l^'Lac) with the 
help of a dichroic waveplate (see the next section). 

It is important that, for better mode matching, beams 
of different wavelengths should be filtered in the trans- 
verse wavevector by different apertures, with the diam- 
eters proportional to the wavelengths. We satisfy this 
condition by separating the wavelengths at the output 
of the MZI, with the help of a dichroic mirror, and us- 



Transformation performed by the dichroic plate. 

The transformation performed by the dichroic plate 
can be described as phase shifts introduced into the cre- 
ation and annihilation operators, different for the two 
wavelengths and two polarizations: 



an - 
bH 



aHe"f^", av 



where 



•JAH 



k°Ad, <t>AV 



ay = ave-'i'^ 
b'y = bye"*'' 



(11) 



(12) 



and similarly for the B beam. Here, k°J^'^g are ordinary 
(extraordinary) wavevectors for beams A, B and d is the 
plate thickness. 

The dichroic plate is a 170/i crystal quartz plate. The 
ordinary and extraordinary wavevectors k'^'^g for quartz 
can be calculated using the Sellmeier equations (see, for 
instance, 01 )• As a result, we get the o-e phase delay for 
Xa equal to (pAH — 't'AV = 4.8547r and for Xb, 4>bh — 
4>BV = 3.7747r, their difference being I.OBtt sa tt. 
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Then the transformation pl|) turns the Hamiltonian 
generating the triplet state l^'mic), 

H' = ihG{a''jjbl + alb''jj) +h.c., (13) 

into 

+h.c. « iftGe*('^-^"+'^«^)(a|j6|, - a[,6|^) + h.c, (14) 

which, up to a phase factor, coincides with the Hamilto- 
nian H generating the singlet state \'^mac)- 
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